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Abstract 

In this note, we discuss the derivation of a formula that has been used 
in the literature in order to compute the number of photons emitted by 
a hot or dense system during a finite time. Our derivation is based on a 
variation of the standard operator-based S-matrix approach. The short- 
comings of this formula are then emphasized, which leads to a negative 
conclusion concerning the possibility of using it to predict transient effects 
for the photon rate. 



1 Introduction 

Electromagnetic radiation (photons and lepton pairs) has long been thought 
to be a good probe of the early stages of heavy ion collisions. Indeed, the 
production rates of these particles increase very rapidly with temperature and 
therefore are dominated by the early times. In addition, photons and leptons 
are interacting very weakly with quarks and gluons, which means that final state 
interactions can in general be neglected. 

The photon and dilcpton production rates have been evaluated completely at 
leading order in the strong coupling a s for a quark-gluon plasma in local thermal 
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equilibrium. These rates have been evaluated at 1-loop 1 in H9QH5I, at 2-loop 
in H3 03 , and the Landau-Pomeranchuk-Migdal ^3 El E] corrections 
have been resummed in El El El El El EH • Up-to-date reviews of the 
situation regarding photon and dilepton production in equilibrium can be found 

in [201 EH E21 123 

The situation is in a far less advanced state when the system is not in local 
thermal equilibrium. Some attempts have been made in order to take into 
account in the calculation the effect of a departure from chemical equilibrium. 
This is done by introducing fugacities in the quark and gluon distributions. 
The fugacity dependent rates have been evaluated at 1-loop in [22], at 2-loop 
in and generalized to the infinite series of diagrams that contribute to 

the LPM effect in These calculations are based on a minimal extension 

of the equilibrium formalism, that merely consists in replacing the equilibrium 
statistical distributions by non-equilibrium distributions. As a consequence, 
they may suffer from the so-called pinch singularities. It can however be shown 
that the corrections due to the "pinch terms" are negligible if the process under 
study is fast compared to the relaxation time |26) . 

In recent years, a new "real-time" approach has been proposed in order 
to compute out-of-equilibrium - namely time-dependent - effects for photon 
production in a dense equilibrated quark-gluon system, which originate in its 
finite life-time |^|[2E|- The result was that there are important transient effects 
that make the yield much larger than what would have been expected by simply 
multiplying the equilibrium rates by the corresponding amount of time. 

This unexpectedly large photon yield, combined with the fact that the energy 
spectrum is not integrable 2 , was the starting point of many discussions regarding 
the validity of this approach [211 021 GH| • In this note, we critically discuss the 
assumptions underlying the formulas used in |27l 1281 To that effect, we 
propose a derivation of the formula giving the photon yield which is based 
on the standard canonical formalism. In particular, we show that it can be 
obtained from an "almost standard" iS-matrix approach, the only departure 
from the standard being that we need to turn on and off the electromagnetic 
interactions at some finite times t{ and tf. 

Our paper is organized as follows. Section 2 is devoted to a derivation of the 
photon yield up to the order e 2 of an expansion in the electromagnetic coupling 
constant. In section 3, we come back to the assumptions made in the derivation 
of the previous section, and discuss their relevance and validity. Section 4 is 
devoted to concluding remarks. 

The loop counting refer to diagrams of the effective theory one obtains after resumming 
the hard thermal loops HI El . 

2 This appears to have changed in more recent iterations 1291 1301 of this work, thanks to a 
subtraction of the terms responsible for this behavior. 
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2 Photon production 



2.1 General framework 

Let us now consider a system of quarks and gluons, and denote by H QCD its 
Hamiltonian (containing all the QCD interactions). We couple the quarks to 
the electromagnetic field in order to study photon emission by this system, 
and denote H e/m , the Hamiltonian of the electromagnetic field, and H ql the 
term of the Hamiltonian that couples the quarks to the photons. The complete 
Hamiltonian is therefore: 

H = Hqcd + H e . m . + H ql . (1) 
At the level of the Lagrangian, this reads: 

£ £qcd £o.m. + t-'qt ! 

C QCD = -\g« v G» v +^(i0 x - gp(x) - m)^(x) , 
r = —-F F^ v 

C qi = -e^{x)A{x)^{x) , (2) 

where G M and G a are respectively the gluon field and field strength, and 
F^ the photon field and field strength, and the quark field (only one flavor is 
considered here), g is the strong coupling constant, and e is the quark electric 
charge. We denote collectively by C lnt the sum of all the interaction terms. 

The number of photons measured in the system at some late time is given 
by the following formula: 

" v £ ^(*><Sw<4 A >)) ■ (3) 

pol. A 

V is the volume of the system and Z = Tr(p(tj)) the partition function. The 
sum runs over the physical polarization states of the photon. p(i;) is the den- 
sity operator that defines the initial statistical ensemble. The "in" states and 
operators of the interaction picture are free, and are defined to coincide with 
those of the Heisenberg picture at the initial time t — U. The "out" states and 
operators are those used in order to perform the measurement. In principle, the 
measurement should take place after the photons have stopped interacting, i.e. 
one should count the photons at a time t — > +oo so that they are asymptoti- 
cally free photons. Here, for the sake of the argument, we are going to define 
the "out" states and fields at some finite time tf. This means that we assume 
that electromagnetic interactions have been turned off before the time tf, for 
this measurement to be meaningful. These "out" states and fields are related 
to the "in" states and fields by means of the "S'-matrix" : 

\a ou t) = S*\a ia ) , 
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a™(k) = S*a£\k)S ) 

aW(k)=Si a W(k)S, (4) 
with an S matrix given in terms of the interaction as 3 : 



S = U(tf,ti) = Vexpi J d xC iat ((f>j n (x)) , 

= U{ti,tf) = Vexpi [VxAntOMaO) • (5) 



V denotes the path-ordering. Expressing the "out" creation and annihilation 
operators in terms of their "in" counterparts, we have: 



dN 

2uj 



4E ^ Uti)StaW(k)a£\k)s) . (6) 



d 3 xd 3 k V V Z 

pol. A 

The next step required in order to bring this expression to an easily calcula- 
ble form is to rewrite the creation and annihilation operators in terms of the 
corresponding fields 4 : 



a Wt (fc) = - ie M*(k) J d 3 xe- ik * dxa A£(x) , 

«L A) (fe) = 4 A) (*0 / d 3 xe ikx d X0 Af n (z) , (7) 



where (k) is the polarization vector for a photon of momentum k and po- 
larization A. The symbol dx is defined as follows: 

A(x ,x) dx B(x , x) = A(x , x) (d Xo B(x , x)) - (d Xo A(x a ,x)) B(x , x) . (8) 

Using the property U (ti , t2)U (t 2 , t 3 ) = U(ti,t 3 ) and the previous relations, we 
can rewrite Eq. © as follows: 



pol. A 

x ^Tr \p(ti )U(t h tf)A& (x , k)A v in (y , -k)U(t f ,ti)} 



1 

Z 

(9) 



where we have performed a Fourier transform on the spatial dependence of the 
photon field: 

A? n (x , k) = J d 3 x e^A^xo, x) . (10) 

3 Note that this is the S-matrix for a system in which the interactions are switched on at 



the time ti and switched off at the the time tj. 

4 The interaction picture field A|^(x) being a free field, these relations indeed give creation 
and annihilation operators that are time independent. One can therefore choose arbitrarily 



the value of xq when evaluating Eq. ■ 
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In Eq. (O, one still has the freedom to chose at will the times xq and yo- This 
freedom can be exploited in order to write: 



2uj 



dN 



cPxcPk V 



pol. A 



x-Tr 
Z 



xo=yo=tf 



(11) 



In this formula, the time path C goes from ti to tf along the real axis, and then 
back to t{. Although it could be included in the definition of C, the portion 
of the contour that goes from tf to +00 and back to tf would not contribute 
to the expression (because U(tf, +oo)J7(+oo, t/) = 1). In Eq. (JTTJ, the time 
derivatives act only on the right (and should be taken before the times Xq 
and yo are set equal to tf). The subscript (— ) (resp. (+)) indicates that the 
corresponding field lives on the lower (resp. upper) branch of the time-path; 
they are necessary in order to ensure that the two fields remain correctly ordered 
when acted upon by the path ordering operator. 



2.2 Order e° 

At this point, one can perform an expansion in the electromagnetic coupling 
constant, while conserving strong interactions to all orders. This is motivated 
by the very different magnitude of the electromagnetic and the strong coupling 
constants. The term of order in the electric charge e is obtained by replacing 
the full interaction Lagrangian £- nrt by the QCD interactions only (which we 



denote £? t °), which leads to: 



2w 



dN 



d 3 xd 3 k 



i J2 ei X) *(k)el x Hk)(d Xo +iu>)(d yo -iu>) 



pol. A 



x(A^-\x ,k)^ + \y ,-k)) , (12) 

\ / x =y =t s 



where we denote: 



Tr 



(A) 



p(U)VAe^c c 



QCD 



(13) 



Tr \p(U)] 

the ensemble average of an operator A at order zero in the electromagnetic 
coupling (but to all orders in the strong coupling constant). Since this object is 
evaluated at order in the electromagnetic coupling constant, the correlator that 
appears in eq. H12[> is nothing but the free path-ordered photon propagator (the 
exponential containing the strong interactions drops out because the photons 
do not couple to quarks at this order 5 ). Therefore, we have in the Feynman 

5 Quarks and gluons can only enter in disconnected vacuum- vacuum diagrams that are zero 
when the time integrations are carried out on both branches of the closed time path C. 
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gauge: 



A^-\x ,k)A^ + \y Q ,-k)) = 



= -<r 



v_ 

2^ 



(1 + ni {uj))e- w{xa - y ^ + n. ( {Lo)e^^ x °- yo) 



(14) 



where we have used explicitly the fact that xq is always posterior to yo on the 
time- path. We denote by n 7 (w) the statistical distribution of photons present in 
the initial ensemble described by p(U). The prefactor V comes from translation 
invariance: the correlator is proportional to a momentum conservation delta 
function (27r) 3 <5(fc — fe) = V. Applying the operator (d Xo + iuo)(d Va — iu) on the 
previous correlation function, we obtain trivially: 



dN 



d 3 xd 3 k 



J2 e^*(k)e^»(k) 



pol. A 



n 7 (u;) = 2n y (uj) 



(15) 



Naturally, this result was expected: at order 0(e°), the number of photons in 
the system is given by the initial photon distribution function multiplied by the 
number of physical degrees of polarization. This is zero if we assume that the 
system does not contain any photons initially. 



2.3 Order e 2 

Expanding cxp — ie f ipfttp to order e 1 produces a 3-photon correlator, which 
vanishes by Furry's theorem. Therefore, the next non zero contribution can 
occur only at the order e 2 . Expanding the exponential containing the electro- 
magnetic interactions up to this order, we obtain 6 : 



2uo 



dN 



d 3 xd 3 k 



e 
~2V 



pol. A 

xU^-\x ,k)A^ + \y ,-k)L 2 ) 



(16) 



with the shorthand: 



A A [ d 1 
du dv / — - 

J (2tt 



d 3 p d 3 q 



) 3 (2tt) 3 



Ai n {uo,p)A( n (v Q ,q)Ji ntP (uo,-p)J inta (vo, -q) . 

(17) 

Jin,p — V'inTpV'in is the electromagnetic vector current. The correlator that 
appears in this formula can now be expanded using Wick's theorem. Naturally, 
the fcrmionic currents can only be paired among themselves. One of the pairings 
corresponds to a term that contains a disconnected vacuum- vacuum factor: such 

6 Implicitly, this is the correct e 2 expansion only if the initial density matrix does not 
depend on the coupling constant e, i.e. if there are no electromagnetic interactions in the 
system at the initial time. 
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a vacuum- vacuum diagram is zero when evaluated on the closed time path C. 
We are left with: 

(A^-\x ,k)^ + \y ,-k)L 2 ) 
d 3 p d 3 q 



duo dvo 



(Jin,p(U0, -p)Jin,*(V0, -<?)) 



(2tt) 3 (2tt) 3 

x{{A^-\xo,k)Ar n (u ,p)){A» n i+) (y ,-k)A° n (vo,q)) 
+<^ ( - ) ( 2 : ,fe)^ 1 ( t ;o, ( 7))(A^ (+) (yo,-fe)Af n ( Mo ,p))}.(18) 

The 2-photon correlator is: 

(A^(x ,k)AUyo,k')) = -(2ir) 3 S(k + k')g^G(x ,y ;k) , (19) 
where we denote 



G(x ,y ;k) = l.[(e c (x Q -yo)+n^u;))e- i ^ Xo - y ^ 

+ (9 c (y -x a ) + n,(u))e-^y°- x « ] } 



(20) 



with 9 c (xq — yo) the generalization of the step function to the time path C. This 
formula is a generalization of eq. (|14|l to the case where the times xq and yo can 
lie anywhere on the time path. The current-current correlator can be written 
in terms of the photon polarization tensor as follows: 

(Jin, P {u , -p)Jin,a{v , -q)) = (2n) 3 S(p + q)U pa (u , v ; -p) ■ (21) 

Acting on eq. Ijl8|l with the operator {d XQ + iu>)(d yn — iuj), we obtain easily: 

dN 



2lo 



d 3 xd 3 k 



= ~E 4 A) *(fc)4 A) (fc) J du dv e 



iuj(vo—u ) 



pol. A 



x[n 7 (u/) + 9 c {u Q - tj)][n 7 (uj) + 6 c (t+ - v Q )}W (u ,v ; k) (22) 

where the superscript + or — on the time tf indicates on which branch of the 
contour C the corresponding time must be considered. 

At this point, we can break down the contour integrations in domains where 
the 8 functions have a constant value. Using the property: 

W^ + (uo,vo;k) + nt l/ _(u ,v ;k) = U t ^_(u ,v ;k) +Ut 1/ + (uo,v ;k) (23) 

which is valid even out-of-equilibrium, we obtain: 



2w 



dN 



d 3 xd 3 k 



e' 
' 2 



J2 4 A) *(fc)4 A) (fc) / du dv < 



,iuj(v — u ) 



pol. A x, 

x[n 7 Hn^_K,« ;fe) - (l + n y (Lu))Tl^ + (uo,vo;k)] .(24) 
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Performing the sum over the photon polarizations, and combining the order e° 
and the order e 2 results, we have the following photon distribution at time tf. 

dN e 2 it 

2uj l3—nT = 2n ^) + "o" / du * dv v e^ v °-^ 
d A xd A k 2 J 

ti 

x[n 7 (w)n£ + _(n ,u ;fc) - (1 + n 7 (w))n^_ + (uo,wo;fc)] 

+0(e 4 ) . (25) 

This formula can be applied to several particular cases, that we discuss in the 
following sections. 



3 Photons, quarks and gluons in thermal equi- 
librium 

A first situation to consider is that of an initial ensemble which corresponds 
to quarks, gluons, and photons in thermal equilibrium at a given temperature 
T. This means that the initial distributions of these particles are given by 
the Fermi-Dirac and Bose-Einstein functions at temperature T. Moreover, the 
initial density operator is p{U) = exp(— H/T) where H is the full Hamiltonian 
of the system, including all the interactions (strong as well as electromagnetic) . 
In this case 7 , the photon polarization tensor obeys the KMS identity, which 
reads: 

n 7 (w)n^l(ito,Vo;fe) = (l + n 7 (w))n^(tt ,uo;fc) ■ (26) 

This relation implies that the e 2 term in eq. (|25|l vanishes. One could in fact 
check that all the higher order terms (in e 2 ) in the photon distribution would 
also vanish thanks to the KMS identity. Naturally, this result is expected: if 
the initial ensemble is an equilibrium ensemble, the populations of particles do 
not change over time. It is also natural that in the calculation this property 
appears as a consequence of KMS: indeed, KMS is the manifestation of thermal 
equilibrium at the level of the Green's functions. In practice, the situation 
considered in this section is only of academic interest, because in applications 
such as heavy ion collisions the photons are not in equilibrium with the strongly 
interacting particles. 



4 Photon-free initial state 

More interesting is the situation of a system which does not contain any real 
photons initially, i.e. for which n 7 (w) = 0. In this case, the photon population 

7 When p(ti) = exp(—H/T), there is an explicit dependence of the initial density matrix on 
the coupling constant e. It is known that this dependence is taken care of without changing 
any formula, simply by appending a vertical branch to the time path |33l . going from t; to 
U —i/T. Moreover, it has been shown that the contribution of this vertical branch can be taken 
into account if one enforces the KMS condition at each intermediate step of the calculation 

|53l 1331 1351 . 
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at time tf is given by: 



2 "d^Pk = ~T f du ° dv ° e " ( " " uo)n ^+("o^o; fc ) • ( 27 ) 



If the initial density matrix describing the distribution of quarks and gluons is 

such that h( u o> vq] k) is invariant under time translation, we can introduce 

the Fourier transform of the photon polarization tensor: 

+00 

n^ + (u ,v ;k) = J — e iE ^~^U^ + (E,k) , (28) 

— OO 

and we can write: 

+00 

— OO 

If one takes the derivative with respect to t/, one obtains the number of photons 
produced per unit time and per unit phase-space at the time tf. 

+00 

2 "dMW = ~ e 2^ njt_ + (^,fc). (30) 



This formula is equivalent to the formula obtained by Boyanovsky et al. How- 
ever, we have derived it here within the framework of an S'-matrix formulation. 
This was not the standard S'-matrix approach though, as some extra hypothesis 
and some extensions have been used. They are discussed in the next section. 



5 Discussion 

A first consequence of eq. ()30|l is that it gives back the usual formula for the 
photon production rate at equilibrium if one takes to infinity the time tf at 
which the measurement is performed (which amounts to turn off adiabatically 
the electromagnetic interactions only at asymptotic times) . Indeed, one has the 
following property: 

lim sin((E-o; ) ft / -t i )) =7r ^_ a;); ^ 



which implies: 



E 



dN e 2 

2w S^ = -T n --+^' < 32 > 
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i.e. we recover in this limit the usual relation between the photon production 
rate and the on-shell photon polarization tensor. In particular, there is no 
contribution to the rate at zeroth order in a s . 

At this point, the main question is whether the finite tf generalization of this 
formula makes sense as a photon production rate, as invoked in - In order 
to discuss this possibility, it is useful to recall and discuss here the hypothesis 
that have been used in order to derive eq. I|30|l . 

• One may be tempted to interpret the number operator al nt a ou t we have 
defined in the third of eqs. Q as the number operator at the time tf for 
photons still interacting with the system, but there is no warranty that 
this definition of the number of photons agrees with the number of photons 
as measured in a detector, precisely because they are not asymptotically 
free states. 

The only possibility to argue safely that this operator indeed counts ob- 
servable photons is to assume that the system does not undergo electro- 
magnetic interactions after the time tf. Another way to state this is to 
say that the object S which appears in eq. I@J is the standard S-matrix 
connecting free "out" and "in" states only if there are no interactions after 

*/■ 

In any case, it is clearly unphysical to keep tf finite: either we are trying 
to measure non asymptotically free photons, or we have to turn off the 
interactions at a finite time t / . 

• A similar problem arises at the initial time. The derivation we have used 
for eq. (|3U[) assumed that there is no dependence on e in the initial den- 
sity operator p{t i ) 1 which is possible only if there are no electromagnetic 
interactions in the initial state. In particular, imposing rt 7 = in the ini- 
tial state is also forbidding electromagnetic interactions before ti (because 
a system of interacting quarks and gluons that undergo electromagnetic 
interactions will necessarily contain photons as well). 

• It is also known that eq. l|3U|l is plagued by very serious pathologies that 
appear as ultraviolet divergences. Firstly, the r.h.s. of eq. (|3U|l turns 
out to be infinite at any fixed photon energy uj due to some unphysical 
vacuum contributions, i.e. processes where a photon is produced without 
any particle in the initial state. Secondly, the remaining terms, even if 
they give a finite photon production rate, lead to an energy dependence of 
this rate which is too hard for being integrable: one would conclude based 
on this formula that the total energy radiated as photons per unit time 
by a finite volume is infinite, which clearly violates energy conservation. 

It was claimed in that the vacuum terms could be discarded simply 
by subtracting to the r.h.s. of eq. (|3U|) the same formula evaluated in 
the vacuum. This indeed has the desired effect, but is a totally ad hoc 
prescription because nowhere in the derivation of the formula appears this 
subtraction term. Or said differently, since eq. (|3 CJ|> is a direct consequence 
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of the definition of eq. © , whatever is wrong with the final formula signals 
a problem either with this definition or with the model. 

Similarly, the authors of [2HI suggested that the divergence that appears 
in the total radiated energy can be subtracted by multiplying the creation 
and annihilation operators used in the definition of the number of pho- 
tons by some wave function renormalization constants. However, no such 
constants appear in the derivation: the operators a ou t , a^ ut can be related 
to their "in" counterparts directly by means of the S'-matrix. 

6 Conclusions 

In this paper, we have obtained a new simple derivation of a formula recently 
proposed in order to calculate the number of photons produced by a system 
during a finite time. The purpose of this derivation is to be sufficiently trans- 
parent in order to exhibit all the underlying hypothesis. In particular, in order 
to obtain this formula, one would have to impose that the electromagnetic in- 
teractions be turned off before the initial time ti and after some final time tf. 
This is clearly unphysical and not surprisingly leads to serious pathologies. Our 
conclusion is that eq. I|30l) does not make sense when keeping a finite tf — ti 
time interval. The transient effects associated to a finite lifetime of the hot and 
dense system, which could enhance the photon production as suggested in |29|. 
cannot be evaluated that way. 
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